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Abstract 



We construct a microscopic model of thermally excited superdeformed states that 
describes both the barrier penetration mechanism, leading to the decay-out transitions to 
, normal deformed states, and the rotational damping causing fragmentation of rotational 

E2 transitions. We describe the barrier penetration by means of a tunneling path in 
the two-dimensional deformation energy surface, which is calculated with the cranked 
Nilsson-Strutinsky model. The individual excited superdeformed states and associated E2 
transition strengths are calculated by the shell model diagonalization of the many-particle 
many-hole excitations interacting with the delta-type residual two-body force. The effect 
of the decay-out on the excited superdeformed states are discussed in detail for 152 Dy, 
143 Eu and 192 Hg. The model predicts that the decay-out brings about a characteristic 
| decrease in the effective number of excited superdeformed rotational bands. 

o 

O ■ 1 Introduction 

O 

Quasi-continuum gamma-ray spectra observed in heavy ion fusion reactions have attracted 
much attention recently as they carry information on structure of rapidly rotating and ther- 
mally excited nuclei. Indeed studies of the ridge and valley structures and the fluctuations in 
the double coincident spectra revealed occurrence of the rotational damping. Namely the col- 
lective rotation of a deformed nucleus becomes a damped motion when the thermal excitation 
energy is provided to the nucleus [Q, §]. This is in contrast to the rotational band structures 
^ ■ known for the levels near the yrast line. The rotational damping is intimately related to 

a basic feature of highly excited compound states that the wave functions of excited levels 
become complex mixture of many-particle many-hole or many quasiparticle configurations 
due to the residual two-body nuclear interaction [||, ||, ||, |j| . 

The quasi-continuum gamma-rays are also produced in the reactions forming superde- 
formed (SD) nuclei || ^, |8|, |9|, |IOfl . However, their properties are much more complex than in 



normal deformed (ND) nuclei. Recent experiments [1C] reveal that the continuum gamma- 
rays in superdeformed nuclei contain more than one components. This complexity arises from 
a characteristic feature that the superdeformed and the normal deformed states coexist in 
the same region of spin and excitation energy. It is known that the observed SD rotational 
bands keep their identity remarkably well even though the SD bands are embedded in a sea 
of compound levels having normal deformation. This is because the SD and ND states are 
separated by a potential barrier in the deformation space. On the other hand, observed SD 
rotational bands terminate suddenly at spin value around 10 — 30/i by decaying to ND states. 
The sudden decay-out of the SD rotational bands are interpreted as a barrier penetration 
phenomenon [11, 12, 14, [l5|, |?], [T(J. As thermally excited superdeformed states relevant for 
the quasi-continuum gamma-ray spectra are concerned, the barrier penetration is expected to 
be more effective. Therefore, it is important to incorporate not only the rotational damping 
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effects (or the complex configuration mixing) but also the barrier penetration in order to 
describe the thermally excited superdeformed states. 

It is interesting to study the thermally excited SD bands because we can learn about 
the collective rotational motion and the collective motion in the shape degrees of freedom 
at the same time. The rotational damping phenomenon is related to the study of quantum 
chaos and is crucial to investigate the very existence of rotation bands in nuclei at thermally 
excited states. The decay-out of the SD bands as a barrier penetration problem tells us 
information on how nuclear shape evolves as a function of the excitation energy and the 
angular momentum. The potential energy surface in deformation coordinates is often used in 
the nuclear structure problems, but the reliability of the energy surface has been tested only 
near the minimum points in most cases. The barrier penetration problem gives a rare chance 
to explore the portions of the energy surface far from the minimum points, i.e. it reflects 
the effect of the large amplitude shape dynamics extending from one minimum to the other 
through a barrier region. The significance of studying the thermally excited SD bands is 
similar to that of spontaneous fissions, but the present problem is unique in the sense that it 
gives possibility to examine such shape dynamics under the influences of the thermal motion 
and collective rotational motion. 

In the present paper, we attempt to construct a microscopic nuclear structure model 
including both the rotational damping and decay-out effects. Theoretical treatment of the 
barrier penetration mechanism relevant to the decay-out of the SD bands was formulated in 
a consistent way firstly by Vigezzi et al. [11]. The validity of the theory have been checked by 
more firm theoretical considerations |Tq] , and also by analysis of the experimental data ]p^] . 
The calculation based on this theory combined with a microscopic potential energy surface 
of the cranked Nilsson-Strutinsky type and the collective mass parameter of the pair hopping 
model gives good account of the sudden decay-out of the SD bands in A 150 nuclei [jD| . 
On the other hand, a microscopic model of the rotational damping in thermally excited 
rotating nuclei has been formulated also on the basis of the cranked Nilsson-Strutinsky mean- 
field ||, H]. This model treats not only the mean- field but also the shell model configuration 
mixing by incorporating the many-particle many-hole excitations built on the cranked Nilsson 
potential and by performing a shell model diagonalization of the residual effective interactions. 
Quantitative success of the model has been demonstrated recently for normal deformed nuclei 
in the rare-earth region |||, || and also in A ~ 110 deformed nuclei [17]. It was also applied 
to the rotational damping in the superdeformed nuclei ]l8| , 18, [9|. However, so far the model 
does not take into account decay-out caused by the barrier penetration. In this paper, we 
extend this approach by combining with the barrier penetration model [11, [l3|] . Previously 
Monte-Carlo statistical simulation models that combine the two effects were proposed [jl4| , |?| . 
However, these models are designed in phenomenological ways so that they can be used 
to analyze the experimental data with parameter fitting or with use of deduced parameter 
values. It is important for further quantitative study to construct a theoretical model on 
the basis of the microscopic description of the rotational damping and the decay-out. The 
purpose of the present paper is to provide with such a microscopic model. 

In the next section we describe details of the model formulation. After recapitulating 
briefly the cranked shell model description of the thermally excited superdeformed states, 
we discuss in detail the barrier penetration mechanism relevant for the thermally excited 



superdeformed states. Here we extend the barrier penetration model of Refs. fll| , 12, |13; [ 
developed for the yrast SD bands in order to describe the decay-out of the thermally excited 
SD states. In §3, we present main results and discuss their implication on the properties of 
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the thermally excited SD states and associated gamma-ray transitions. The calculations are 
done for typical superdeformed nuclei, 152 Dy, 143 Eu and 192 Hg in the A ~ 150 and 190 mass 
regions. In §4 we summarize the points of this paper. 



2 The Model 

2.1 Configuration mixing and E2 transitions among SD states 

Excited superdeformed states lying highly above the yrast states may consist of many-particle 
many-hole configurations built on the superdeformed mean-field, but they hardly have pure 
mean-field configurations due to the two kinds of mixing effects. One is the mixing among SD 
states caused by the residual two-body interaction. This brings about the rotational damping, 
which leads to fragmentation of the rotational E2 transitions among SD states. The other 
is the mixing between SD and ND states caused by the barrier penetration. This leads to 
decay-out transitions to ND states. In this subsection we briefly recapitulate a microscopic 
formulation of the first kind of configuration mixing among the SD states [Q, 18, 



We start with the cranked Nilsson single-particle Hamiltonian h w = /^Nilsson — ojJ x whose 
quadrupole and hexadecapole deformation parameters are determined by the minimization 
of the potential energy surface calculated by the Strutinsky method. The cranked Nilsson 
Hamiltonian h u defines the single-particle orbits associated with the rotating superdeformed 
potential as a function of the rotational frequency u>. Many-body shell model basis config- 
urations are then defined in terms of the cranked Nilsson single-particle orbits. We employ 
all many-particle many-hole configurations but with a truncation with an upper limit on the 
excitation energy. Changing the variable from the rotational frequency to the rotational spin 
/, where / means the expectation value of J x (angular momentum along the rotation axis), 
the shell model basis states are constructed for each spin and parity, as usually prescribed by 
the cranking model. This shell model basis for the superdeformed states is denoted 

In the shell model Hamiltonian describing the mixing among {\fj,(In))}, we include the 
volume-type delta force v (1, 2) = v T 5{x\ — £2) as the residual two-body interaction. The force 
strength v T is taken the same as Refs. |l8|, |l9| . Diagonalizing the shell model Hamiltonian, 
we obtain a set of eigen solutions for each spin / and parity it. We denote their energy E a (In), 
and state vectors \a(Iir)) = ^2^X^{I-k) \/j,(Iit)). The B(E2) strength of stretched rotational 
E2 transitions from an eigenstate a at spin / to states (3 at / — 2 are calculated by using the 
obtained wavefunctions. It is convenient to normalize the B(E2) strength so that the sum 
over final states for a fixed initial state becomes the unity. We denote the normalized strength 
S a l, /3I -2, which approximately corresponds to the branching ratio of the E2 transitions from 



the state |a(/7r)). See Refs. |18| , 19, |4J for details of the cranked shell model calculation. 



The calculation of \a(Iir)) , E a (Iir) and 5 , Q / j/ g/_2 in the present paper is the same as those in 



Refs. 1 18, [19]. 



2.2 Barrier penetration and decay-out transitions to ND states 

Because of the barrier penetration effect, the superdeformed states \a(In)) obtained in the 
cranked shell model will couple further to normally deformed (ND) compound states which 



lie energetically near E a {IiT). We follow the theory proposed by Vigezzi et al. [11] to describe 
this coupling between SD and ND states. 

It is assumed that ND compound states are distributed randomly around a SD state 
|a(/7r)), and that the coupling matrix elements between |a(/7r)) and the ND states have an 
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average size of v. The Vigezzi model describes this coupling in terms of a quantum tunneling 
taking place through the potential energy barrier in the deformation space. This leads to the 
relation between v and the tunneling width r t as, 

r t = — - = 2nv 2 Pn , (1) 

where p n is the level density of the ND states and D n = l/p n is the mean level spacing. 

Due to the coupling, some components of normal deformed states are mixed in the su- 
perdeformed state \a(Iir)). Since the normal deformed states have also electromagnetic tran- 
sition probabilities, the SD states thus mixed have possibility to decay not only with the 
rotational E2 gamma-rays feeding other SD states but also with the gamma-rays feeding 
to normal deformed states. The latter is the origin of the decay-out transitions. For the 
electromagnetic transitions associated with ND states, we consider the statistical El and the 
rotational E2 transitions. A key quantity that reflects the effect of decay-out is the branch- 
ing ratio of the decay-out transitions, N out ; it depends on the transition probabilities of the 
competing electromagnetic transitions, i.e. the E2 transition width r s among superdeformed 
states, and the El and E2 transition width r n for the normal deformed states, and also de- 
pends on the amplitudes of the normal deformed states mixed in the superdeformed states. 
The mixing amplitudes may fluctuate strongly depending on the relative positions between 
the SD and ND states. In addition statistical fluctuation in coupling matrix elements is 
caused by the compound nature of the ND states. We can evaluate the average value of iV ou t 
assuming random matrix model description of the ND states and the Gaussian orthogonal 
ensemble (GOE) [jll|. The resultant average value (N out ) of the decay-out branching ratio is 
a function of the tunneling width It or the coupling matrix element v, the decay width r n of 
the ND transitions, the decay width r g of the SD transitions, and the average level spacing 
D n of the ND states. More specifically, (N out ) is a function of the two ratios r t /D n and r s /r n 
in the model of Ref. [11|, and we have calculated its functional form by numerically diago- 
nalizing GOE matrices of dimension four hundreds. Recently theoretical models to calculate 



(-^Vout) have been proposed [£l], 22, 16]. Approximate treatment of Ref. [21] gives a simple 
result, (iVout) = A/ (A + r s ), but it was found that the approximation is not valid in the 
realistic situations of decay-out. In Ref. |]l^] the statistical theory using the supersymmetry 
technique has been rigorously solved in general conditions, and a simple two-level model has 



been investigated in [22], which is appropriate in the weak coupling limit. We examined the 



result of the Vigezzi model [11] by comparing with that of [16|, and we found that both 



theories give the same results within the relevant range of parameters r s ,r n ,r t and D n for 
the decay-out of SD bands in both the A ~ 150 and 190 regions. 

If the tunneling width i~t, the level spacing D n of the ND states, and the electromag- 
netic decay rates r s and r n are given, one can calculate the average decay-out probability 
{N ut(E, -0) for any specified values of the spin and the excitation energy of the superdeformed 
states. It is then possible to include the effect of the decay-out on the excited superdeformed 
states \a(In)) calculated by the shell model configuration mixing. Now the state |a(/7r)) with 
energy E a has an average branching ratio 1 — (N on t(E a , I)) for the rotational E2 transitions 
feeding to other superdeformed states at I — 2, and the ratio (N out (E a , I)) for the decay-out 
transitions. Thus E2 transition probabilities 5 , Q / j/ g/_2 calculated by the cranked shell model 
(see §2.1) are renormalized as 



S a I,pi-2 =U - (N ou t{E a , I)} ^ S a i : f3i-2- 
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Hereafter we denote the average decay-out probability iV ou t without the symbol of average 
for simplicity of notation. 

In the next section, we analyze properties of rotational E2 gamma-rays from the excited 
SD states in terms of this renormalized E2 transition probabilities. In the rest of this section 
let us describe more details of the model constituents. 



2.3 Collective potential and mass tensor 

In order to specify the tunneling width r t , we assume that the barrier penetration takes place 
as the nucleus changes its shape from the super deformation to the normal deformation. We 
consider the two quadrupole deformation parameters (e2>7) of the cranked Nilsson model in 
the Lund convention as the dynamical variables relevant for the deformation change. The 
collective Hamiltonian for the deformation variables is given by 

1 2 

^=2 E m iMii 92) 4i4j + V(q 1 ,q 2 ) (3) 
ij=l 

where the collective coordinates q± and q 2 are defined by 

qi = e 2 cos (7 + 30°), q 2 = e 2 sin (7 + 30°). (4) 

The collective potential energy V(<7i,<72) is calculated for each signature and parity by 
means of the cranked Nilsson-Strutinsky method with the monopole pairing residual inter- 
action as a function of the deformation parameters (€2,7), or equivalently (qi,q 2 )- Both 
the static and dynamic pairing correlations are taken into account in terms of the cranked 
Hartree-Bogoliubov method and the random phase approximation (RPA) as has been done 
in Ref. |23|. Namely the response function technique is used with small imaginary part of 
frequency, 5 = Im(u>) = 200 keV/K, which is slightly larger value but saves computation time 
considerably. The cut-off parameter for calculating the effective RPA pairing gap |23|, 
which is used in Eq. (S), is taken as u; cu t = 400 keV/h. Note that exchange contributions are 
excluded in both the RPA correlation energy and RPA pairing gap (<5-Erpaj ^rpa in l^ij] )- 
The angular momentum is conserved in the tunneling problem so that the potential energy 
surface should be calculated with fixed spin values. Spin-interpolation is used for this purpose. 
The hexadecapole deformation (54 in the Nilsson potential) is determined so as to minimize 
the potential energy for each point of (£2,7) deformation parameters and the spin values. 
The grid of calculation for the deformation parameters (qi,q 2 , £4) are —0.18 < q\ < 0.60, 
-0.12 < q 2 < 0.42 and -0.10 < e 4 < 0.14 with an interval of 0.06. The grid for rotational 
frequency, which is used for the spin- interpolation, is determined by Auj = AI/J str with 
AI = 3, where J s t r is Strutinsky-smoothed moment of inertia calculated for each set of de- 
formation parameters (e2 ; 7>e4) ; and the maximum frequency is extended to give calculated 
spin value up to 60 h. The Is and I 2 parameters of the Nilsson potential is taken from Ref. pl| . 
The remaining parameter for the potential energy calculation is the monopole pairing force 



strength, which are determined by the smoothed gap method [26| for each set of deformation 
parameters with the model-space energy cut-off of \.2T\uq below and above the Fermi surface 
(hujQ = 41.0/^4 1//3 MeV is the oscillator frequency). Note, however, that in Ref. the 
minimization of energy with respect to the 64 deformation has been done only approximately, 
namely the value determined without pairing correlations has been used, and the strength of 
pairing correlations has been determined by the standard smoothed pairing gap A = 12/y/A 
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MeV. We have found that the full minimization against €4 including the static and dynamic 
pairing correlations leads considerably smaller height of barrier compared with the previous 



calculations [13|, which makes the tunneling probability too large. We, therefore, employed 
smaller pairing force strength determined by A = 10/y/A MeV in the present calculations. 
This value is chosen to obtain barrier height which gives similar tunneling probability to that 
of the previous calculations. The pairing force strength thus determined is smaller by about 
6% for nuclei in both the Dy and Hg region. 

As for the collective mass tenser m^gi, (72), we use the one based on the pair hopping 



model 27, Pq ], where the semiclassical (Fermi-gas) estimate for counting the number of 
pair-level crossing has been extended to the two dimensional case of the volume-conserving 
ellipsoidal deformation [12, The mass tensor is given explicitly asQ 



where An^ and Ap U> are the effective RPA pairing gap p3, 24] for neutron and proton 



in. 



i(eff) 



2Al 



A^\q,I) 



+ 



A^\qJ) 



(hop) 
2 m ij 



(<?), 



(5) 



calculated at each deformation (q) and spin /, and mf* (q) is written analytically as 



m ( ^\q) 



(hop) 



1 da da 1 
+ 



db db 1 dc dc 
+ 



a 2 dqi dqj b 2 dqi dqj c 2 dqi dqj 



with 



in 



(hop) 



h 2 G Q A 2 
A 2 27' 



(6) 



(7) 



Here the three axis lengths, a, b, and c, for the volume-conserving ellipsoidal shape can be 



(hop) r- 

mo Oij, i.e. 



111 



(hop) . 




IS 



specified by the deformation parameter (€2,7), and mf^ op \q = 0) 
the hopping mass at the spherical shape. The empirical values for the constant parameters, 
Go = 25/ A MeV and Aq = 12/\/ r A, are employed p8| . The same tunneling model consisting 
of the collective potential and the mass tensor explained above has been used for the problem 



of decay of high-X isomers, see Ref. [ 30 1 for details. 



2.4 Tunneling decay width 

If the decay-out of the yrast SD band is concerned, the tunneling width of this state can be 
evaluated semiclassically as 



Tluj s 
~2^ 



(1 + 6XP2S)- 1 , 



(8) 



where T = (1 + exp25) 1 is the tunneling transmission coefficient and S is the (imaginary 
time) action integral along a classical tunneling path, 



S(E) 



path 



dsyj2M (V(q(s))-E), 



(9) 



where s is a parameter along the classical tunneling path, and is defined according to 

M ds 2 = ^2mij(q)dqidqj. (10) 



1 There is a misprint of missing factor (3/2) in Eqs. (22) and (26) of Ref. J2!|. It is corrected in Eq. 
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Here Mq is a mass unit and we take Mq = ml hop '. See §2.5 for determination of the path. 



When we apply to the yrast SD band, the energy E of the classical tunneling path is chosen 
as the zero-point energy of the shape motion around the SD minimum. The factor hu) s /27r 
in Eq. @ is the knocking probability where uj s is the vibrational frequency of the collective 
motion along the tunneling path near the SD potential minimum. Eq. (g) has been used for 
description of the decay-out of SD rotational bands [11, 12, [R^, fj], [l5|| . 



The above expression, however, may not be directly applied to the tunneling of the excited 
superdeformed states with high thermal excitation energy since it assumes that the internal 
nucleon configuration is kept adiabatically along the tunneling process, considering only the 
collective shape dynamics. If the highly excited SD states are concerned, the collective motion 
associated with the tunneling will be coupled with the complex internal excitations, which is a 
basic feature of highly excited compound states caused by the residual two-body interaction. 



Following Bj0rnholm-Lynn [31], which discusses the tunneling of fission isomer states, we 
consider that collective vibrational excited states form doorway states of tunneling, and the 
strength of the vibrational states are spread over many compound states. If we assume 
that the coupling is sufficiently strong and the strength is distributed uniformly, the average 
strength of doorway states in a compound energy level is evaluated as ~ D s /Huj s . Here 
D s = l/p s is the average level spacing of the superdeformed states at a given excitation 
energy (p s is the level density of SD states). The knocking probability for the tunneling 
of the compound level is then reduced as ~ (huj s /2ir)(D s /huj s ) = D s /2ir in average. The 
tunneling width in this situation is given by 

r t = ^(l + ex P 2S)- 1 . (11) 

Here the action S of the classical tunneling path is given by Eq. (|9|) whereas the energy E of 
the path is chosen as the energy of the excited SD state under consideration. The knocking 



probability D s /2tt has been used also in the previous simulation model of Ref. [14, pTOf] . 

An important feature of Eq. ( [TT] ) is that it can describe the limit of statistical mixing when 
applied to the excited states located above the barrier height. Since the SD and ND states 
above the barrier are expected to mix completely with the statistical weight {p n vs. p s ) and to 
form the compound states, an average strength of a SD state |SD) in the compound states \i) 
is given by | (i|SD) | 2 p s /(p Q + p s ) = D tot /D s , where D tot denotes the average spacing of the 
statistically mixed compound states. On the other hand, if we assume S = above the barrier, 
the tunneling width given by Eq. ( JTT| ) reads It = D s /Att. Using this expression, the average 
strength of the superdeformed state can be evaluated as | (i\ST>) | 2 « D tot /r t « <iitD tot /D s , 
which agrees with the value for the statistical mixing except a numerical factor. If the factor 
1/2-k in Eq. (0) were removed, better agreement with the statistical limit would be achieved. 



Since the argument of Ref. [31] for the derivation of Eq. (11) concerns an accuracy of an order 
of magnitude, such possibility may also be allowed. As shown later, however, the difference 
in the numerical factor has only minor effect. It is also interesting to see that an expression 
of the tunneling width similar to Eq. (|TT|) is derived for the tunneling of a chaotic system 



with two degrees of freedom on the basis of the semiclassical theory [32|. There the tunneling 
width is proportional to the average level spacing of the initial states and the tunneling 
transmission coefficients, as in Eq. (|ll|), but it also has an additional prefactor related to a 
second order dynamical correction p2| . Eq. ([ll]) does not include such a dynamical factor, 
which is, however, not known for many-body fermion systems. 

We evaluate the tunneling width i~t of the excited SD states by means of Eq. (|ll|) since it 
describes the average value of the tunneling width r t , assuming the uniform distribution of 
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doorway strength. However, if the strength distribution of vibrational states is not uniform 
and concentrated to some excited SD states, the tunneling width for those states may be 
larger than that given by Eq. (|Tl]), and some other states may have smaller tunneling width. 
Although it is not easy to evaluate such fluctuations, we can evaluate an upper limit of r t . 
It is given by Eq. (||) and calculated at the excitation energy under consideration since this 
equation corresponds to a special situation where a excited SD state contains large amount of 
the vibrational strength. We may also use Eq. (@) for the excited SD states in the following 
analysis in order to estimate an upper limit of r t . We adopt Eq. (||) also to describe the 
decay-out of yrast SD states. It is important to note that the two evaluations, Eqs.(|8|) and 
(|TT|), give different results. Since Eq. ( |TT| ) contains the knocking probability D s /2ir that 
decreases steeply with increasing excitation energy, it gives significantly smaller value of the 
tunneling width than Eq. (||) at high excitation energy. As far as the the SD states near the 
yrast states are concerned, on the other hand, the two expression gives more or less similar 
value since the level density D s near the yrast line may be order of several hundred keV, 
which is approximately the same as the value of Hlu s (See Table Q). In the following section, 
we will discuss the difference by using numerical results. 



We remark that it is also possible [33] for the tunneling width i~t to be reduced from the 
average value given by Eq. (|iT| ) , if the strength of doorway states is not uniformly distributed 
but concentrated in between the excited SD states and so the coupling between them is 
minimal. However, we do not evaluate such reduction since a precise modeling of the doorway 
states is required and it is beyond the scope of the present paper. 

2.5 Tunneling path and classical action 

Examples of the calculated Nilsson-Strutinsky potential energy surface V(q\, 92) are shown in 
Fig. |l|. The local minimum at large quadrupole deformation 62 ~ 0.4 — 0.6, 7 ~ corresponds 
to the superdeformed state. The minimum with smaller quadrupole deformation corresponds 
to normal deformation. The potential barrier and the saddle-point separating the two classes 
of states are also seen in the figure. The quantum tunneling takes place through the barrier 
from the excited superdeformed states situated at the superdeformed minimum to the normal 
deformation. 

The tunneling path and the associated classical action integral is calculated in the fol- 
lowing way. For a given energy E, the deformation space (91,92) is separated into classically 
allowed and forbidden regions, whose borders are given by energy contours with energy E. 
There are two allowed regions, one surrounding the SD minimum and the other corresponding 
to normal deformation. The tunneling path is the one that starts at a point on the border of 
the allowed region at the SD side, travels through the barrier region, and reaches a point on 
the border of the other allowed region at the ND side (See Fig.fi]). 

The tunneling path should minimize the action S(E), Eq. (|9|), under a constraint of fixed 
energy E. The minimal action path obeys a classical equation of motion 

v - d 2 qj x ^ drriij dqj dq k 1 x dm jk dqj dq k dV _ 

dr* 8q k dr dr 2 % dr dr % ' 1 ' 

where r = it is the imaginary time and it is related to the path parameter s, Eq. (|i"0|), through 



dr = I Mq 

ds V2{V{q)-EY { ' 
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Namely, the tunneling motion corresponds to the classical motion in the inverted potential 
energy, — (V{q)—E). Solving the classical equation of motion with an initial condition starting 
at a point on the border at the SD side, a classical trajectory is obtained. The equation is of 
the second order and the initial condition is given by qi = en and dqi/dr = at r = 0, where 
the point P^ = (01,02) is on the energy contour at the SD side. The zero initial velocity 
(in imaginary time) is the consequence of the energy conservation. An arbitrary classical 
trajectory does not necessarily reach the energy contour at the ND side, but if we choose a 
proper initial point p(°\ the classical path reaches the ND energy contour. This path, called 
the escape path, gives the minimum action. We search such a solution by changing the initial 
point P( ) of the classical path on the energy contour. This method to find the minimal action 
path is parallel to that of the Schmid's method [34], which deals with the path starting at 
the potential minimum. Since in this method the initial point P(°) is placed at a potential 
minimum and the force dV/dqi = 0, an infinitesimal initial velocity should be imposed in 
some direction, which should be chosen properly to obtain the escape path reaching to the 
contour of the ND side. The search of direction of the infinitesimal initial velocity in the 
Schmid's method just corresponds to the search of the initial point on the contour in our 
case of the finite excitation energy. The numerical calculation is relatively easy in our case 
since it is a one dimensional search in the two dimensional collective coordinate. Depending 
on the shape of the barrier and the mass tensor, more than one solutions may exist as local 
minima of the action. In such cases, we adopt the solution which gives the smallest action. 

An example of the obtained tunneling path is shown in Fig. for the excitation energy 
E ex = 1.5 MeV measured from the SD potential minimum. For energy below the saddle-point, 
the tunneling path and the action S is obtained by the procedure just described. Above the 
saddle-point energy we set 5 = 0. When the excitation energy under consideration becomes 
lower than the ND yrast, there should be no tunneling, i.e., It = 0. In the following, we 
define the barrier energy Vb as a excitation energy of the saddle-point, i.e. Vb = Saddle 
since the action vanishes when the excitation energy exceeds the barrier energy, i.e. S = if 
E >V\). Note, on the other hand, that the least action path does not necessarily go through 
the saddle point, see e.g. Ref. @ and also Fig. g. 

We take into account of the zero-point oscillation of the shape dynamics for the description 
of the tunneling. This can be achieved just by assuming that the yrast superdeformed state 
has an zero-point oscillation energy huj s /2 above the SD minimum (with energy i?g D ) of 
the potential energy surface V(q\,q2). Here uj s is the frequency of the oscillation at the SD 
potential minimum given by a small amplitude approximation along the tunneling path at 
the energy E ~ 22g D (see Table [l]). Namely, the tunneling path and the action S(E) for the 
yrast superdeformed states are calculated with the energy E = £ S Dyr = For the 

excited states having the excitation energy U (measured from the SD yrast state), we calculate 
the tunneling path with the energy E = E^ d + ?iuj s /2 j i-U. In the following, we use U to denote 
the physical excitation energy measured from the SD yrast energy -EsDyr = Eg D + hu! s /2 while 
we use another notation E ex = E—E^y) to represent the energy relative to the potential energy 
minimum Eg D . 

In Ref. [13] a brute force method has been used in order to obtain the minimum action 
path, i.e. the path is divided into finite segments and the positions of all segments are varied 
in such a way to realize the minimum action. In Ref. the Schmid's method has been 
used to obtain the minimum action path at E ex = 0, along which the action is calculated 
with using E ex = %oj s /2 in Eq. ([]). This treatment can be justified if E is much smaller than 
the energy of the saddle point. In the present paper, on the other hand, we solve the least 



9 




ecos(y+30°) ecos(y+30°) 



Figure 1: The calculated potential energy surfaces V{q\,q2) for 152 Dy (/ = 40 + ), 143 Eu 
(I = 61/2+) and 192 Hg (I = 30 + ) shown in the top, the left bottom, and the right bottom 
panels, respectively. Interval of energy contours is 0.2 MeV and 7 = is indicated by an 
oblique line. The SD minima are located around €2 = \J q\ + q\ = 0.45 ~ 0.6, 7 f» 0. In the 
top panel, the energy contour with E = i?g D + 1.5 MeV is shown with thick solid curves, and 
the least action path at this energy with a thick dashed curve. The solid circles indicate local 
minima of the ND side whereas the saddle-point is also shown with a diamond. 
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Figure 2: The tunneling action S(E, I) plotted as a function of spin for different excitation 
energies for 152 Dy (a, n) = (0, +), 143 Eu (1/2, +) and 192 Hg (0, +). The signature and parity 
(a, 7r) for each nuclei is taken the same as the yrast SD configuration, as in Fig. |. The thin 
solid curve with the largest value of action corresponds to the zero excitation energy E ex = 0.0 
(the energy of the SD potential minimum) while the others are for E ex = 0.5, 1.0, 1.5, • • • MeV. 
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0.50 
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0.93 


0.52 
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Table 1: The calculated barrier parameters, hui^, Vb — -EsDyr and huj s , for several spin values. 
The signature and parity is taken the same as in Fig. |. Here I = I — a in the first column 
is a combination of the spin / and the signature a. The unit for the barrier parameters is 
MeV. 

action path explicitly, starting from the finite energy contour as is explained above, since the 
excitation energy ranges up to the saddle-point energy. 

In Fig. H we plot the action value of the classical tunneling path as a function of angu- 
lar momentum for different values of excitation energies for the same parity and signature 
quantum numbers as those in Fig. |l|. The thick solid line corresponds approximately to the 
action for the yrast SD band. (This is because the zero-point energy is about 0.5 MeV.) The 
increase in the action with spin is caused by increasing barrier height, Vb — -^SDyr; compared 
to the SD yrast and by increase in the collective mass due to the quenching of pairing. The 
action S(E) decreases with excitation energy at a fixed spin since length of the tunneling 
path becomes shorter. It becomes zero as the excitation energy U reaches the barrier height. 

The decrease in S(E) with excitation energy is almost linear as seen in Fig. |2[ If the 
potential barrier were one-dimensional and had a inverted quadratic shape, the action integral 
would have a form S(E) = tt (Vb — E)/hu>b where Vb is the barrier energy and Wb is an 
oscillation frequency characterizing the barrier. Using the linear energy-dependence and the 
formula of the harmonic approximation, we can characterize the excitation energy dependence 
of the microscopically calculated action with the two parameters Vb and Wb- In practice, the 



11 




Figure 3: The energies of the ND and SD yrast states, -EwDyr and -EsDyrj and the saddle- 
point energy Vb as a function of angular momentum for 152 Dy (left), 143 Eu (middle) and 
192 Hg (right). The signature and parity is taken the same as in Figs, [l] and ||. 



barrier energy V\> is evaluated as the saddle-point energy, and the barrier frequency Wb is 
evaluated as ftwb = 7r (Vt> — E)/S(E) at the energy E = -EsDyr of the SD yrast state (or 
at the energy of the ND yrast state if the latter is higher). Table |l] shows the calculated 
values of the barrier frequency ftwb and the relative barrier height Vt> — -E-SDyr for several spin 
values, together with the zero-point oscillation frequency uj s around the SD minimum. The 
barrier height has rather strong spin dependence which is common to other calculations of 
potential energy surfaces 36]. Note that the barrier frequency Wb also shows sizable spin 
dependence. This is related to increase of the collective mass caused by decrease of pairing 
correlation at higher spins. 

The energies of the yrast superdeformed and normal deformed states, -EsDyr and -ENDyrj 
and the height of barrier (the saddle-point energy) , Vb , are plotted in Fig. || as a function of 
spin for the parity and signature quantum numbers of the yrast SD state. The correction for 
the zero-point oscillation energy is included in the ND yrast energy I?NDyr in the same way 
as for the superdeformed yrast states: The zero-point energy is taken the same as that of the 
superdeformed minimum. 

Here difference of potential barrier in Fig. |3| among these nuclei should be emphasized. 
The barriers in 152 Dy and 143 Eu increase monotonically (more precisely, the barrier height of 
143 Eu almost saturates at I ~ 40ft, which is clearly seen from the action in Fig. [2] for E cx <; 1.5 
MeV), while that in 192 Hg increases at low-spins but starts to gradually decrease at I ~ 30ft. 
This fact is seen directly in Table |], and the value of actions precisely reflects the effect of 
this feature of barrier height in Fig. ^. Namely, the action of 192 Hg does not increase linearly 
as function of spin, while that of 152 Dy does, especially at higher excitation energies; the 
action even changes to decrease around I ~ 25ft as the barrier height does. This feature that 
the barrier height saturates and changes to decrease at I ~ 20 — 30ft is commonly observed 
for nuclei in the A ~ 190 region, which can be seen in other calculations 36]. We will see 
in the next section (see §||) that the decay-out properties in 192 Hg are dramatically different 
from those of 152 Dy and 143 Eu. 



2.6 Level density 

The average level spacing D n and D s of the ND and SD states are calculated as follows. Using 
the same cranked Nilsson single-particle orbits as employed in the shell model diagonalization, 
we construct excited many-particle many-hole configurations for each spin and parity. The 
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Nilsson quadrupole deformation parameters are fixed to those associated with the normal 
deformed (or superdeformed) minimum. The density of these excited configurations are 
calculated as a function of the excitation energy by explicitly counting the configurations 
in energy bins. For normal deformed states the resultant level density as a function of the 
excitation energy fluctuate but nicely follow in average the Fermi gas formula appropriate for 
the cranking model for fixed spin and parity |37 



PFG (U) = ^a-^U-Uxp2^U. (14) 
This allows us to use the Fermi gas expression to represent the average level density p n 



of the normal deformed states. The formula used in Ref. 13] was the one like Eq. (14) 
but with factor two being multiplied. This was because the states with different parity mixes 
through the statistical El transitions. However, the level density in our model is used through 
Eq. ([!]), where the mixing is caused by the tunneling process which conserves the spin and 
parity. Therefore, we believe that Eq. ( |i~4|) should be used. 

As for the level density parameter a, we use the one that is adjusted to fit the microscopic 
calculation of the level density. Since spin dependence of the calculated level density is found 
to be weak, we adopt a single value of a for all spin values by neglecting the spin dependence. 
The fitted value is a n = A/9.43, A/9. 27 and A/10.15 for normal deformed states in 152 Dy, 
143 Eu and 192 Hg, respectively. The values of the level density parameter is different from 
the traditional values a = A/8 or A/10. This is because of the shell effect in the density 
of the single-particle orbits near the Fermi surface. These values are close to a recent semi- 
empirical estimate a n (emp.)= A/8. 57, A/9. 46 and A/11.32 for 152 Dy, 143 Eu and 192 Hg based 
on the neutron resonance data |3q| . The average spacing is given by D n = 1/ pFG,n(U f , I) by 
using the fitted Fermi gas formula Pfg,ii(U' , I), where the excitation energy U' is the energy 
measured from the yrast ND energy. We take the lowest ND state among the four set of 
the parity and signature quantum numbers as the yrast configuration to define the yrast ND 
energy, for which the zero-point energy is included. For the other signature than that of the 
yrast configuration, we make an interpolation to define the yrast ND energy. 

We adopt a similar procedure for D s of the superdeformed states. Since the excitation 
energy dependence is slightly different from Eq. (|i"4| ) due to the strong shell effect, we make a 
fit at the excitation energy E cx = 2.5 MeV which is relevant for the quasi-continuum spectra. 
The obtained level density parameter is a s = A/18.9, A/18.5 and A/13.0 for 152 Dy, 143 Eu 
and 192 Hg 



2.7 Electromagnetic decay widths 

We evaluate the electromagnetic decay width in a way similar to Ref. |jl3; ]. For the SD states, 
only collective rotational E2 transitions are considered, and the decay width is given by 

r s = Tf 2 D = 3.0 x 10~ 4 Q 2 E^ (15) 

in the unit of MeV where the quadrupole moment Q (in fm 2 ) is calculated by the liquid 
drop model for the deformation parameters of the SD minimum. The transition energy E~ is 
evaluated for each spin and parity as the energy difference between the SD potential minima 
at / and 1 — 2. We use the same value for the excited SD states by neglecting excitation 
energy dependence which are expected to be small. As for the effect of rotational damping, 
we note that the rotational damping widths for the thermally excited SD bands are less than 
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100 keV [19], while the E2 transition energy is typically 500 — 1000 keV. Therefore it is also 
expected to affect the electromagnetic width only slightly and so its effect is neglected. 

For the transitions associated with ND states, the collective rotational E2 and the sta- 
tistical El transitions are taken into account. Namely the decay width is given by r n = 
r E2 D + r Ei°- Here r E2° is calculated in the same way as in the SD case. In ND case, how- 
ever, there often exist more than one local minima with different quadrupole deformations. 
Therefore, we take average of the E2 width over the local minima. The weight factor for 
averaging is taken proportional to the level density associated with each local minimum eval- 
uated with the fitted Fermi gas formula pFG,n- This averaging avoids possible discontinuity of 
-Tg 2 D associated with crossing of different yrast configurations. In addition, even if the lowest 
minimum is spherical i^P remains finite. The averaging, however, has little influence on the 
result except in the case in which the yrast minimum turns out spherical. Dependence of 
-Tg 2 D on excitation energy is weak as far as the yrast state is well deformed. The El decay 
width is calculated by using the GDR strength function /gdr> 



t-iND 
1 El 



Cei 



u> 



PfgAU') 



(16) 



where U' is the energy of the SD state relative to the ND yrast, and Cei is the El hindrance 
factor. Actually we use the analytic expression for which is a good approximation 

at excitation energy U' under consideration, 



7-1ND 
1 El 



Cei 2.3 x 10~ n NZA 1/3 



a 



5/2 



(17) 



where the Thomas-Reiche-Kuhn sum rule value is used for the strength function /gdr and 
the GDR parameters, -Tgdr = 4 MeV and £gdr = 80A -1 / 3 , are employed. We use the El 
attenuation factor Cei = 0.15 taken from Ref. [4C]. The El width T 
fifth power of the temperature, T = y/U'/a with increasing the excitation energy U' . 



increases with the 



3 Results 

3.1 Decay-out of yrast superdeformed bands 

We have performed calculations for yrast SD bands in both the A ~ 150 and A ~ 190 regions, 
in which the calculational procedure has been improved over the previous calculations [ jDj] , 
as is outlined in the previous sections. Here we use Eq. (||) for the tunneling width i~t, 
appropriate in describing the tunneling of the yrast SD bands. The only definite observable 
for the decay-out of the yrast SD band is the 7-ray intensity, of the transitions I — » I— 2 

inside the band. This quantity is calculated by using the average decay-out probability iV ou t 
discussed above as 

J 7 (J) = J] (l-iV out (C/ = 0,/ / )). (18) 

r>i+2 

We show in Fig. ||, the calculated 7-ray intensities by solid lines for selected yrast SD bands 
in 152 Dy, 143 Eu and 192 Hg in comparison with the data; the intensity data are taken from j|l]] 
for 152 Dy, from p2|, El for 143 Eu (see also [44j for the most possible spin-assignment), and 



from [45| for 192 Hg. In contrast to the previous calculations [13|, where schematic rotational 
spectra were used for both the SD and ND yrast bands, all the quantities entering in the 
calculations are obtained theoretically; namely, there is no adjustable parameters in the 
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Figure 4: The 7-ray intensity for decay-out of the yrast SD band in 152 Dy (left), 143 Eu 
(middle) and 192 Hg (right). The solid, dashed and dot-dashed lines denote the results without 
any adjustment, with modification by C p , with modification by C mass (both are explained in 
the text), respectively. The result with using the "back-shifted level density" with E^s = 1.0 
MeV, is also included as a dotted line for 192 Hg. 



present calculations. As a result, it is not so easy to reproduce the observed intensity. The 
pattern of the 7-ray intensity is characterized in two ways; one is the decay-out spin I ut, 
that is defined as a spin value at which the 7-ray intensity is reduced to half, i 7 (i ou t) = 1/2, 
and the other is a sharpness of the intensity-drop. The sharpness is well reproduced in all 
three SD bands as was in the previous calculations [13], especially if the decay-out spin is 
adjusted in some way, see the following discussions. The problem lies in the calculated decay- 
out spin, which is too large in 192 Hg, i^t — ~ H^-> while it is rather well reproduced 
in 152 Dy, l^ut ~~ ^out P ' > ~ 3ft, and in 143 Eu, m —2ft. In order to set up a reliable model for 
the thermally excited superdeformed states, it is very important to reproduce the decay-out 
spins. Therefore, let us examine possible origins of the above discrepancies. 

As is discussed in the previous sections (see §2.2), the decay-out probability, iV ou t, and 
so the 7-ray intensity, i 7 , is determined by four quantities; the tunneling width, it, the 
electromagnetic decay widths for the SD and ND bands, r s and -T n , and the mean level 
distance of ND states, D n (or equivalently the level density, p n ). Table Q shows the values 
of these quantities responsible for the decay-out, evaluated at XJ = 0, I i^f an d I^nt^ ■ 
As for the electromagnetic transition rates, r s is that of the stretched E2 and r n is of the 
statistical El in the decaying spin region. They are well determined by the equilibrium 
shape and the excitation energy of the SD band from the ND yrast band, namely by the 
minima of the calculated potential energy surface. On the other hand, it is calculated by 
the collective mass tensor and the potential energy barrier, i.e. the behaviour of the saddle 
point region. Our potential energy surface gives correct equilibrium shape, which is deduced 
from the measured life time of the SD band, and the calculated barrier also seems to be 



consistent with other calculations, see for example Refs. [35, 36]. Therefore we do not try 



any adjustment for quantities related to the potential energy surface in the present work, 
although there may exist a room for some improvement if we use a different mean-field from 
the Nilsson potential or a different type or strength for the residual pairing interactions. 
Remaining quantities in question are the level density and the collective mass parameters. 
We try to adjust these quantities in order to obtain agreement of the decay-out spin as in 
the following. 



15 







152 


Dv 




143 


Eu 




!92 Hg 




1=28 


1=32 


1=37/2 


1=41/2 


1=12 


1=22 




5.9 x 10" 


-5 


4.7 x 10-° 


3.7 x 10" 


-5 


1.1 x 10" 


5 


5.9 x 10 


-5 


2.8 x 10" 


-6 


D n 


1.1 x 10" 


-5 


1.2 x 10" 5 


5.9 x 10- 


-4 


1.1 x 10- 


3 


5.7 x 10" 


-6 


2.2 x 10" 


-5 


r s 


1.8 x 10" 


-8 


2.9 x 10~ 8 


1.5 x 10- 


-9 


3.2 x 10- 


9 


5.8 x 10- 


-11 


2.9 x 10" 


-9 




5.7 x 10" 


-9 


5.6 x 10~ 9 


1.3 x 10- 


-9 


1.0 x 10- 


9 


4.3 x 10 


-9 


3.7 x 10- 


-9 


"sDyr 


4.68 




4.61 


2.57 




2.25 




4.24 




3.53 




S 


4.0 




5.2 


4.3 




4.9 




4.0 




5.4 




'S'(Cmass) 


5.1 










3.7 




8.2 









Table 2: The calculated values of i~t, D n , r s , r n and the action S at U = 0, I « I^ ut and 
-^out P ^ > responsible for the decay-out of the yrast SD band. The energy U^^of the SD yrast 
states measured from the ND yrast is also listed. The last line lists the value of the action S 
at / ~ -^out P ' ) calculated with the collective mass renormalization C mass = 1.4, 0.7 and 3.0 for 
152 Dy, 143 Eu and 192 Hg, respectively. The unit for r t ,D n ,r s ,r n and C/^ Dvr is MeV. 



We remark that some additional effects on the level density are not included in the present 
calculations. For example, the pairing correlation that may be present at relatively low spins 
may reduce the level density by a factor of up to ~ 10 2 while the effects of collective vibrations 
such as beta and gamma modes may increase the level density. To represent these effects and 
some other origins of ambiguity, we shall allow a renormalization factor of the level density 
C p that multiply the level density p n of normal deformed states as p n — ► C p p n . Since we 
explicitly deal with the superdeformed states in the shell model description, we do not apply 
such a renormalization for the superdeformed states to avoid an inconsistency. The results 
fitting the decay-out spin by using the factor C p are shown by the dashed lines in Fig. the 
values of which are C p = 0.1, 15 and 2 x 10~ 4 , for 152 Dy, 143 Eu and 192 Hg, respectively. It 
should be mentioned that in the calculation of Ref. the calculated 7-ray intensity was 
mistaken by two spin units, i.e. it was — 2); this is the main reason why the C p factor 
is about an order of magnitude smaller in the present calculation than in |i~3T| . The resultant 
C p factor in 152 Dy and 143 Eu may be acceptable. But it is too small for 192 Hg since the ND 
level spacing at the excitation energy U' ~ 4 MeV becomes D Q ~ 30 — 100 keV, which is 
unphysically large. The fluctuation analysis of the quasi-continuum decay-out transitions in 
192 Hg [46 1 suggests that the ND level density p n is smaller than the Fermi gas estimate due 
to the pairing correlation, and that the reduction of p n is the order of 10 _1 — 10~ 2 . As for 
the level density, the back-shifted Fermi gas formula, namely U is replaced by (U — Ebs) m 
Eq. (14), is often used to include the effect of the pairing correlation. In the decay of the Hg 
region, the decay-out spin is small and this effect may be more important than in the case of 
the Dy region. Therefore, the results with using Ebs = 1.0 MeV is also included as a dotted 
line for 192 Hg in Fig. ||. As seen in the figure, it helps to shift the calculated decay-spin lower 
by about 3 K but the amount is not enough. 

As an another origin of the discrepancy, we consider a possibility that the collective mass 
could be affected by some additional effects. The essential contents of the pair hopping 
model is the number of single-particle level crossings and the pairing gap. We have estimated 
the number of crossings by using a lowest order semiclassical (Fermi-gas) approximation in 
Eq. @, and detailed properties of individual nucleus cannot be reflected. However, if one 
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calculate the number of crossing microscopically or take into account the shell effects in 
the single-particle orbits, this quantity could be largely modified. This may influence the 
collective mass Tnij(qi,q2)- To represent such effects, we introduce a renormalization factor 
C m ass that modifies as rriij(qi,q2) — > C mass rriij(qi, 52)- It should be noticed that the scaling 
factor C mass affects the calculated action in two ways: one is a simple scaling of action by a 
factor VCmass as is clear from Eqs. (|9|) and (|l0|), the other is a change through the variation 
of the zero-point energy hui s /2 in the action integral. Both effects make the action larger 
(smaller) if C mass > 1 (< 1). The results fitting the decay-out spin by using the factor C mass 
are shown by the dot-dashed lines in Fig. the values of which are C mass = 1.4, 0.7 and 3.0, 
for 152 Dy, 1 Eu and 192 Hg, respectively. As in the case of the introduction of C p , the factor 
seems large for Hg. 

From the above discussion, it is clear that we have large discrepancy for the decay-out 
spin value in 192 Hg. According to the results of more systematic calculation [47] it is one 



of the worst cases in the A ~ 190 nuclei. In this systematic calculation the empirical level 
density parameters [38] are used so that the results for the intensity data are slightly different 



from those shown in this paper. Also, the back-shifted level density with £bs = 1-0 MeV is 
used for comparison with the data for nuclei in the A ~ 190 region. The results of Ref. [^] 
show that the basic characteristic features of intensity pattern are reproduced in both the 
A ~ 150 and 190 regions; especially the common rapid decrease of transitions at lower spins. 
The observed decay-out spin is about 25h on average in the A ~ 150 region, while it is 
about 10?i in the A ?a 190 region. On average calculated decay-out spins are comparable 
to the observed ones in the A ~ 150 region, but those in the A ~ 190 region are a little 
larger than the observed ones systematically even though the back-shift in the level density 
formula is employed. The calculation reproduces qualitatively the empirical observation that 
the average decay-out spins in the A ~ 190 nuclei are systematically smaller than those in the 
A ~ 150 nuclei, but the amount is not enough. This trend is mainly due to the facts that the 
calculated action as a function of spin is somewhat different between nuclei in the two regions 
(see Fig. ||), which reflects the difference in the potential energy surfaces, and that the level 
density parameter, a n /A, is systematically smaller in the A ~ 190 region. We have to admit 
that the calculation overestimates systematically the decay-out spins in the A ~ 190 region, 
which suggests that the calculated action values are too small in these nuclei. It seems that 
some further mechanisms which discriminate SD bands in these two regions are necessary. 
Presently we do not have any clue to solve this problem. 

3.2 Decay-out of excited superdeformed states 
3.2.1 Tunneling width 

Let us first discuss the tunneling width It which is the most important quantity in describing 



the decay-out of excited superdeformed states. As discussed in §2.4, we evaluate the tunneling 



width with use of Eq. (11), which is expected to give an average behaviour of It- In the 
analysis below, we shall also show results obtained by using Eq. (^), which provides an upper 
limit of It, in order to illustrate the difference in the two evaluations of It- The use of 
Eq. ( |TT|) is referred to as the case (a) in the following, while the evaluation of the upper limit 
with Eq. @ is as the case (b). 

In Fig. H we show dependence of the calculated tunneling width on the spin and the 
excitation energy. The upper three panels show results in the case (a), calculated for the 
three nuclei, 152 Dy(left), 143 Eu(middle) and 192 Hg(right). The calculated tunneling width It 
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Figure 5: The tunneling width r t for 152 Dy (left), 143 Eu (middle) and 192 Hg (right) as a 
function of spin for different excitation energies. The line styles are defined in the same man- 
ner as in Fig. g. The upper three panels show the results in the case (a) using r t = (D s /2tt)T 
whereas the lower panels show those in the case (b) calculated with i~t = {huj s /2ir)T, The 
signature and parity is the same as in Fig. pi 



increases exponentially with the excitation energy until the energy reaches the barrier height. 
Note that the tunneling transmission coefficient T = 1/(1 + exp25) increases exponentially 
with excitation energy since the action integral S decreases almost linearly with increasing 
the excitation energy (cf. §2J5 and Fig. ^). On the other hand, the knocking probability 
proportional to the level spacing D s exhibits an exponential decrease, but this decrease is 
weaker than the increase in the transmission coefficient. A systematic spin-dependence is also 
seen in It for the energy region below the barrier height. Namely, the tunneling width for 
the yrast SD states (corresponding to the thick curve in Fig. ||) decreases exponentially with 
increasing spin except at high spins (I ^ 30) in 192 Hg. The spin-dependence originates mainly 
from the tunneling transmission coefficient since the spin-dependence in the level density 
of SD states is weak. As the excitation energy exceeds the barrier height, the calculated 
tunneling width turns to decrease. This is because the transmission coefficient saturates to 
a constant value ~ 1/2 and the tunneling width r t ~ D s /4ir then decreases with excitation 
energy as it is governed by the level spacing D s of the excited SD states. In this regime, the 
tunneling width depends mostly on the excitation energy measured from the SD yrast, but 
it is insensitive to the spin as the level density p s of SD states is. 

The lower panels show the tunneling width in the case (b) where the knocking frequency 
is given by the collective vibrational frequency u s . Since u; s depends little on the excitation 
energy and spin, the energy dependence of the tunneling width is solely governed by the 
exponential increase of the transmission coefficient below the barrier, and it saturates at 
i~t ~ htu s /4:Tr above the barrier and becomes spin independent. 
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Figure 6: The average decay-out probability N out (E,I) plotted as a function of spin, calcu- 
lated for 152 Dy ((q, n) = (0, +)) with C p = 0.1. The thick solid, dashed and dot-dashed lines 
indicate the results in the case (a) using It = (D s /2ir)T for excitation energies E ex = 1.0, 
2.0 and 3.5 MeV, respectively. Thin lines represent the results in the case (b) employing 
r t = (Huj s /2it)T. 

Comparing the two cases quantitatively, it is seen that the tunneling width in the case 
(a) is significantly smaller than that in the case (b) by more than the order of 10 2 for the 
excitation energy E ex = 2 ~ 3 MeV and the difference increases further as E ex reaches near 
and above the barrier. (Note, however, that the value of it in the case (b) for the energy 
above the barrier height is not meaningful since Eq. (||) does not give correct description of 
the statistical mixing that is expected above the barrier.) We see also that the difference of 
it between the two cases is not significant at the low excitations near the yrast SD states 
(E ex ^~> 1 MeV). This is because D s in this energy region is of the same order as htv s ~ 1.0 
MeV. 



3.2.2 Decay-out properties 

Figure |6| shows the calculated average decay-out branching ratio iV ou t in 152 Dy. The value 
of iV ou t is plotted as a function of spin I for different excitation energies E ex = 1.0,2.0 and 
3.5 MeV measured from the SD potential minimum -Eg D . A steep increase of -/V" ou t with 
decreasing I around iV ou t ~ 1/2 indicates a sharp decay-out, and this feature is commonly 
seen for all excitation energies. As the excitation energy increases, the decay-out occurs at 
higher spins. This is because the tunneling probability increases. In Fig. [6L we also show 
iVout calculated with use of Eq. (||) (the case (b)) with thin lines. The difference between 
the two cases (a) and (b) is significant for high excitation energy E ex ^ 2 MeV because of the 
difference in it becomes large as discussed above. On the other hand, the two evaluations 
give similar results for E cx 55 1 MeV. 

The decay-out branching ratio N out modifies the E2 transition probability of exited su- 



perdeformed states as described in §2.4. The energy levels E a (I) of SD states and the E2 
transition probabilities S a j^j-2 calculated for 152 Dy are plotted in Fig. [7[ In this figure 
shown are such strong E2 transitions that satisfy S a j : pj^2 > 1 /V% by connecting the energy 
levels with solid lines. Compared with the previous calculation which does not include the 
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Figure 7: The energy levels of the yrast and excited SD states, shown with horizontal bars, 
calculated for 152 Dy ((a,7r) = (0,+)) by means of the cranked shell model diagonalization. 
A rigid body rotational energy E vig ^ = I {I + l)/i7 r i g id is subtracted. Solid lines connecting 
the energy levels indicates strong E2 transitions satisfying S a i^i-2 > 1/V%- 

decay-out effect (Fig. 6 in Ref. ]18[), strong E2 transitions disappear at low spin and high 
excitation energy due to the decay-out to normal deformed states. The sequence of strong E2 
transitions associated with the yrast SD states correspond to the yrast SD rotational band, 
and the figure shows that the yrast SD band decays out around I = 28 (See also Fig. ||). 
There exist several rotational bands above the yrast band up to the excitation energy U ^ 2 
MeV. The decay-out spin of these bands is higher than that of the yrast bands since the 
decay-out branching ratio N out increases with excitation energy. 

Since the decay-out occurs at the spin where iV ou t increases sharply from zero and reaches 
to N out ~ 1/2, it is possible to define the decay-out spin I ut(U) for a given excitation energy 
U by means of a condition N out = 1/2. This definition of the decay-out spin is slightly 
different from that used in § |3.1| to describe the decay-out of the yrast SD band (It was 
defined there by ri/'(l ~~ -^out(-O) = 1/2)- However, the actual difference is less than 2h 
for the value of the decay-out spin, and in the following discussion we adopt the condition 
-Wout = 1/2 as the definition of the decay-out of excited superdeformed states. Since N out is 
a function of the spin / and the energy E, the criterion iV ou t = 1/2 gives a boundary in the 
(E,I) plane. We can represent the decay-out boundary also in terms of the decay-out energy 
E OVLt (I) as a function of spin / as well as the decay-out spin I out (U). 

The decay-out energy E out (I) calculated for 152 Dy is shown in Fig. ||. Results in the left 
panel is obtained with no adjustment of the level density nor of the mass parameter. In 
this figure, we show not only the results in the case (a), evaluated by using r t = (D s /2n)T, 
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Figure 8: Left panel: Decay-out energy E out (I) for 152 Dy ((a,7r) = (0,+)) as a function 
of spin i, together with the barrier height Vf, and the ND yrast energy E-^D yr . All these 
quantities measured the SD yrast energy -EsDyr are plotted. The decay-out energy are shown 
in the three cases, (a) it = (D s /2ir)T, (a') it = D S T and (b) it = (hu) s /27r)T. The 
calculation is done without any adjustment of the ND level density and the collective mass. 
Middle panel: The onset energy of rotational damping -Edamp(i) is plotted together with the 
boundary energy for the SD-ND mixing defined by r t /D n = 1 and the decay-out energy E out , 
calculated in the case (a) with C p = 0.1. The three regions separated by the two boundaries 
-E-out(i) and -Edamp(i) are indicated as "pure rotational band" , "damped rotation", and "ND 
dominant compound" (see text). Right panel: Same as the middle except for the adjustment 
of the collective mass with C mass = 1.4. 



Eq. (|Tl"l), which gives the average behaviour of the decay-out, but also two other results in 
the cases (a') and (b). In the case (a') we assumed it = D S T that gives the tunneling width 
different from that in the case (a) by a factor of 2n. Comparing the results in the two cases (a) 
and (a'), we see that the ambiguity in the numerical factor of the tunneling width (difference 
of the factor 2tt) does not cause large difference in the decay-out boundary. In the case (b), 
on the other hand, we evaluate with use of the other expression it = (htv s /2ir)T, Eq. (g), 
which is expected to give an upper limit of i~t. The decay-out energies E out (I) for the two 
cases (a) and (b) differ by several hundred keV at maximum, and the difference is larger than 
the ambiguity of a numerical factor in the cases (a) and (a'). We can conclude from these 
comparisons that the reduction of the knocking probability (the factor ~ D s /lu s in Eq. ( |j~l| ) 
v.s. Eq. (H)) plays important roles for the description of the excited SD states. 

The decay-out boundary for the zero excitation energy U = corresponds to the decay- 
out of the yrast SD band. As the excitation energy increases, the decay-out boundary moves 
to higher spins. In other words, the decay-out of excited superdeformed states takes place 
at higher spin than that of the yrast SD band. The calculation indicates that the decay-out 
spin I ou t(U) increases by about 8h for increasing excitation energy by 1 MeV. Equivalently 
the decay-out energy E ont (I) increases monotonically with increasing spin. The slope of the 
decay-out energy E out (I) is larger than that of the barrier height Vt>(i) since the decay-out 
probability iV ou t depends not only on the tunneling transmission coefficient but also on other 
quantities such as the level density of ND states, which vary strongly with spin. 

The decay-out spin for U = is seen in the figure to be I ou t(U = 0) « 30?i, which is 
slightly larger than the experimental decay-out spin of the yrast SD band. In order to have 
more reliable description, we have discussed in § |3.1| two ways of adjustment by introducing 
the renormalization factor C p of the ND level density and by the renormalization C mass of 
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the collective mass parameter. In the middle (right) panel, we show the decay-out boundary 
for the excited SD states, evaluated in the case (a), with use of the adjusted ND level density 
with C p = 0.1 (the adjusted collective mass parameter C mass = 1.4). Since the renormaliza- 
tion factors C p and C mass are not very different from unity, difference between the left and 
middle (right) panels is not large. The difference between the two adjustments (middle vs. 
right panels) is very small for the lower energy region U ^ 2 MeV since the decay-out spin 
of the yrast SD band is adjusted in both cases. The difference, however, is not negligible 
for higher energy region, especially for the region above the barrier height. The renormaliza- 
tion of the collective mass affects the decay-out property only below the barrier, while the 
renormalization of the ND level density influences both below and above the barrier. 

In order to clarify the decay-out properties, it is useful to recall that the average decay-out 
probability N out is a function of two ratios r t /D n and r s /T n |ll|, |DJ. We found that the 
behavior of the decay-out energy E out is mainly governed by the ratio r t /D n . To show this, 
we plot in the middle and right panels of Fig. || a boundary energy that is defined by the 
condition r t /D n = 1, and compare with the decay-out energy E ont (I). If the ratio r t /D n is 
more than 1, it implies that a significant amount of components of ND states mix into SD 
states. The condition r t /D n = 1 defines a boundary for the strong mixing among SD and ND 
states. It is seen from the figure that the boundary of strong mixing follows rather well with 
the decay-out energy E ont (I)- The ratio of the electromagnetic transition widths r s /r n plays 
a relatively minor role in determining the spin and excitation energy dependence of the decay- 
out, especially for the spin region much higher than the decay-out spin of the yrast SD band. 
For the high spins, the E2 decay width of ND states dominates over the El decay width, and 
the ratio of electromagnetic widths r s /r n is approximated by (Q s /Qn) 2 x {Jn/Js)^, where 
Qs (Qn) and J s (J7n) are the static quadrupole moment and moment of inertia for the SD 
(ND) states. The ratio then depends neither on excitation energy nor spin, and takes a value 
which is a little more than 1 but does not exceed 10 1 . 

The decay-out boundary E on t(I) in the region above the barrier height shows behaviour 
different from that below the barrier. As seen in the middle panel (and also in Figs. ^ for 
143 Eu and 1C for 192 Hg), the boundary rises very sharply above the barrier. In the region 
above the barrier the SD and ND states mix statistically to form compound states as the 
ratio rt/D n reduces to the ratio of the level spacings D s /D n of SD and ND states. (Note that 
the results of the cases (a) or (a') should be adopted for the region above the barrier since the 
case (b) is not applicable in this region.) As spin increases, the ND yrast energy relative to 
the SD yrast becomes higher and the ratio r t /D n ~ D s /D n = p n /p s decreases exponentially 
for a fixed value of U. In other words, superdeformed states become dominant components 
of the statistically mixed compound states. On the other hand, the ratio r s /r n increases 
with spin since the E2 decay probability r s associated with SD states becomes larger due 
to the E^ factor while the ND decay width r n decreases due to the decrease of the relative 
excitation energy. Both two effects, the decrease in r t /D n ~ D s /D n and the increase in 
r s /r n , favor dominance of the SD transitions at very high spins while the ND transitions 
dominates at lower spins. The decay-out spin above the barrier height represents a border 
where the relative dominance of the SD and ND transitions is interchanged. 



3.2.3 Rotational damping effects 

Let us now discuss effects of the rotational damping on the excited superdeformed states. 
As seen in Fig. 0, The strong E2 transitions with S a i^i_2 > 1A/2 gradually disappear in 
the region with excitation energy U greater than about 2 MeV. This disappearance of the 
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strong E2 transitions is caused by the rotational damping. Although the transition from 
the rotational band structures seen near the yrast line to the rotational damping in the 
highly excited states is gradual, let us characterize it by introducing a boundary energy 
for the onset of rotational damping. For this purpose, we calculate the branching number 
^branch = (J2a &ai /3/-2) _1 o HI' EOI that quantifies the fragmentation of E2 decay from a 
SD state a at spin I. The onset energy of damping -Edamp(-0 is then defined as the energy 
where the average branching number calculated as a function of the energy E exceeds 2 for 



a given spin /. The calculation of -Edamp(-0 is the same as in Ref. [1^, 19] . The onset energy 
-^damp(^) of the rotational damping is plotted in the middle and right panels of Fig. |8[ 

Using the decay-out energy E out (I) and the onset energy -Edamp(-0 of rotational damping, 
the excited superdeformed states are divided into three regions. The E2 transitions associated 
with the superdeformed states can exist only in the low energy region below the decay-out 
boundary E out (I) or at the high spin side of this boundary. This region is further divided into 
two by the second boundary -Edamp(-0- Below E^ amp (I), the superdeformed states form the 
rotational band structures, i.e., sequences of levels connected by the strong E2 transitions. 
The superdeformed rotational bands are confined in this region, which is labeled in Fig. || 
with "pure rotational band". On the other hand, concerning the excited SD states which 
are located above -Edamp(-0 and at the high spin side of the decay-out boundary E ovl t(I), 
they decay also by the rotational E2 transitions characteristic to the SD shape, but the the 
rotational E2 transitions are fragmented (i.e. the rotational damping) due to the admixture 
of superdeformed many-particle many-hole configurations caused by the residual interaction. 
This second region is labeled in Fig. [|| with "damped rotation" . Note that the damped E2 
transitions can be found even above the barrier. The third region that is situated at low-spin 
side (and high energy side) of the decay-out boundary E out (I) is labeled with "ND dominant 
compound" in Fig. |8[ The SD states in this region mix strongly with the ND states, and 
the decay of these states are dominated by the El or E2 transitions associated with the ND 
states. 

We may use the onset guideline of whether the excited superdeformed 

states have character of compound states. At excitation energy higher than -Edamp(-0> the 
SD states become strong admixture of different SD np-nh configurations. Thus it may be 
justified to use the average estimate of the tunneling width, given by Eq. ( [il|) (the case (a) 
with the knocking frequency oc D s ), for the energy region sufficiently higher than -Edamp(-0- 
But for the energy region below -Edamp(-O) we can expect possible fluctuation around the 
average estimate. As discussed above, an upper-limit estimate of the fluctuation is given by 
Eq- (|8|) (the case (b)). Note, however, that the ambiguity due to the fluctuation is not large 
in this energy region since the difference in the decay-out boundary E ovl t(I) between the two 
cases is small (see the left panel of Fig. ||). 

3.2.4 143 Eu 

Figure || shows the results for 143 Eu. Basic features are found similar to those in 152 Dy. 
The influence of the adjustment of the ND level density (the middle panel with C p = 15.0) 
and of the collective mass parameter (the right panel with C mass = 0.7) on the decay-out 
boundary is not large, likewise in the case of 152 Dy. There exits also in this nucleus the three 
regions classified in terms of the decay-out and the onset of the rotational damping, as shown 
in the middle and right panels. A remarkable difference from 152 Dy is that the decay-out 
boundary is located at lower spin by about Wh. The low values of decay-out spin is seen not 



only for the yrast SD states U = (as discussed in §3.1) but also for the highly excited SD 
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Figure 9: Same as Fig. | but for 143 Eu ((a,ir) = (1/2,+)) with C p = 15.0 (middle panel) 
and with C mass = 0.7 (right panel) 




Figure 10: Left panel: Calculated decay-out energies E out (I) in the cases (a) and (b) (see 
text) for 192 Hg ((a,ir) = (0,+)) calculated with no adjustments [C p = l,C mass = 1), plotted 
together with the onset energy of rotational damping Ed& mp (I), the boundary energy for 
the SD-ND mixing defined by r^jD n = 1 calculated in the case (a), the barrier height Vb 
and the ND yrast energy E^Dyr- Middle panel: The decay-out energies E out (I) in the case 
(a), calculated with the level density renormalization C p = 2.0 x 10~ 4 . Right panel: The 
decay-out energies E out (I) in the case (a), calculated with the collective mass renormalization 
Cmass = 3.0 and the boundary energy for the SD-ND mixing. 



states including the states above the barrier. Consequently the region of "damped rotation" 
is much larger than in 152 Dy. The calculation predicts the decay-out spin J out ~ 35 — 45h 
for the highly excited SD states above the barrier height. The SD states lying Ji> 35 — 45h 
and U ^ 2 MeV emit characteristic rotational E2 gamma-rays with fragmentation caused 
by the rotational damping. This feature seems consistent, at least qualitatively, with the 
experimental observation of intense quasi-continuum E2 transitions pj. 



3.2.5 192 Hg 



Figure 1C shows the results for 192 Hg. Although there are some common features to those 
in 152 Dy and 143 Eu, large differences are observed. The left panel shows the results without 
any adjustment. It is seen that difference between the case (a) and the case (b) is larger than 
in 152 Dy and 143 Eu. This is because the spin dependence of the tunneling width r t is much 
weaker than in 152 Dy and 143 Eu as seen in Fig. [|, hence the difference in It (the cases (a) 
vs. (b)) causes larger change in the decay-out boundary. Weak spin-dependence is also seen 
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in the barrier height. These features are due to the qualitative difference of the calculated 



potential energy surface between the Dy and Hg nuclei, which is stressed at the end of §2.5. 
It is also noted that the decay-out boundary above the barrier is located at rather low spin 
(I rs 40h). This is caused by the relatively high level density of the super deformed states in 
this nucleus. (The corresponding level density parameter is a=A/13.0 MeV -1 .) An unwanted 
feature in the left panel is that the theoretical prediction of the decay-out spin for the yrast 
SD band is significantly higher than that observed in the experiment by more than 10h as 
already discussed in § [3,1| . 

In §[0], we introduced two ways of adjustment by renormalizing the level density of ND 
states by a factor C p = 2.0 x 10~ 4 and by renormalizing the collective mass parameter by 
C m ass = 3.0. The results of these calculations are shown in the middle and right panels. 
Remarkably large difference between the two calculations is seen for the decay-out boundary 
of the excited SD sates although both are adjusted to reproduce the observed decay-out of 
yrast SD band. The change in the collective mass parameter increases the action value of the 
tunneling path by a factor of 1.9—2.1 and reduces the tunneling transmission coefficient, but 
it does not influence the behavior above the barrier where the tunneling action becomes zero. 
On the other hand, the change in the ND level density influences both below and above the 
barrier. Since the renormalization factor for the adjusted ND level density (C p = 2.0 x 10 -4 ) 
is very small, its effect is drastic. For example, the decay-out boundary above the barrier is 
largely shifted to I « 20ft, which is very different from the decay-out boundary calculated 
with the adjusted collective mass (C mass = 3.0). Note, however, that the modification of the 
ND level density with the factor of 10 -4 may not be realistic, although it is also not certain 
whether the renormalization of the collective mass parameter by a factor of 3 is justified. We 
can say that the decay-out boundary has different spin dependence and stays at considerably 
lower excitation energy in 192 Hg compared to 152 Dy and 143 Eu if a realistic ND level density 
is used, and this is mainly due to the fact that the barrier height in the A ~ 190 nuclei 
saturates or even decrease gradually at 20 — 30ft. 

3.3 Number of rotational bands 

As Fig. |7] illustrates, the number of superdeformed rotational bands existing in a nucleus is 
limited to a finite number since the rotational E2 transitions associated with highly excited 
levels exhibit the rotational damping. In addition, the decay-out transitions are present at 
low spins. Since the rotational band structures are not observed if the decay-out probability 
dominates over the E2 transitions, the decay-out further reduces the number of rotational 
bands especially at low spins where the decay-out plays a role. An effective number of decay 
paths, which approximately gives the number of superdeformed rotational bands, has been 
evaluated experimentally from analysis of the quasi-continuum E2 gamma-rays in superde- 



formed nuclei y, |9|, ^8[. In contrast to the intensity data for the decay-out of yrast SD bands, 
where the tunneling probability in a narrow spin-range is relevant because of the sharp decay- 
out, the number of rotational band is sensitive to the decay-out process in a broad spin-range, 
as is shown in the following, and is effective to study the barrier penetration problem. We, 
therefore, adopt in this subsection the number of rotational bands as a typical observable 
on the excited superdeformed states, and discuss in detail effects of the decay-out on this 
quantity. 

The number of superdeformed rotational bands, iVband> is evaluated as follows. To make 
a distinction between the levels forming rotational bands and those exhibiting the rotational 
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Figure 11: The number of SD rotational bands, A^and) f° r 152 Dy, calculated in the two cases 
(a) and (b) (see text) with C p = 0.1 and plotted as a function of spin /. Also shown are the 
result with C mass = 1.4 in the case (a) and the one without the decay-out effect. 



■>P SaI,/3I-2 



-1 



damping, we consider the branching number of E2 transitions branch = (j2 
defined for each SD state a at spin /. The effect of the decay-out on this quantity is taken 
into account by means of Eq. (||). Namely, the branching number is now given by 
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"-branch (a, I) 



aI,pi-2 



(1 - N out (E a ,I)) Ti b ranch(a, I) 



(19) 

(20) 



by multiplying the factor (1 — N out (E a , I)) 2 to the original branching number Tibranctr 

The 

modified branching number takes a very large value if the decay-out transition is dominant 
(i.e. iV ou t ~ 1). The condition n-branch < 2 gives a definition of the states having strong E2 
transitions forming superdeformed rotational bands. For fixed spin and parity, we count the 
excited superdeformed levels that satisfy n-branch < 2 for consecutive two steps of strongest E2 
transitions [|, 18, [OJ. Summing over four combinations of parity and signature quantum 
numbers, we evaluate the number of superdeformed rotational bands Aband- 

We show in Fig. 11 the number of superdeformed rotational bands calculated for 152 Dy. 
Here we plot -/Vb an d obtained with the renormalization of the ND level density {C p = 0.1) and 
of the collective mass parameter (C maS s = 1.4) that are chosen to reproduce the decay-out spin 
of the yrast SD band (cf. § |3.l[) . The average estimate of the tunneling width r t = (D s /2tt)T 
(the case (a)) is adopted. The result without the decay-out is also plotted for comparison, 



which is the same as that obtained in the previous paper [18]. We find that the effect of decay- 
out reduces significantly the number of bands for the spin region I ^>50h. This behavior is 
easily understood from Figs. and ||. Since the decay-out boundary cuts the region of pure 
rotational bands at I ~ 26 for U = (yrast states) and at / ~ 50 for U ~ -Edamp, the effect 
of the decay-out on -/Vb an d become sizable in the spin region / ~ 26 — 50h. At high spins 
7^50 where E out (I) exceeds -Edamp^), the decay-out has little influence on ./Vband- There 
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Figure 12: Same as Fig. 11 but for 143 Eu. Here iVband is calculated either with the adjusted 
C p = 15.0 or with the mass renormalization C mass = 0.7. It is plotted as a function of 
the average gamma-ray energy E~, and compared with the experimental effective number of 



decay-path [48]. 



is only minor difference between the two ways of renormalization. In the figure we plot also 
iVband obtained with the upper limit estimate of i~t = (huj s /2ir)T (the case (b)). The case (b) 
can be considered to provides a lower limit for the evaluated -/V^and- The difference in iVband 
between the two cases (a) and (b) is up to about by 10 ~ 15. The number of bands is as 
much as 50 — 60 for long spin interval / ~ 45 — 60. It is emphasized that this number is larger 
than those for ND rare-earth nuclei (iVband ~ 30), indicating that the shell effect leading to 



the large value of A^, an d 19] remains even with the presence of the decay-out effect. 
Figure 12 shows the results for 143 Eu. Here the experimental data ]f48l are compared. To 



make a direct comparison, we plot A^and as a function of the E2 gamma-ray energy. The 
transformation from the spin to the gamma-ray energy is done by calculating the average 
gamma-ray energy calculated over all E2 transitions as a function of spin. We show the 
calculated results for both cases (a) and (b). The difference between (a) and (b) and also the 
difference between the two renormalizations are found smaller than that in 152 Dy, and it is 
not more than the size of experimental error bars. The decay-out has little influence on the 
number of bands for spin region / <; 50 (£" 7 1400 keV). Gradual decrease of number of bands 
with decreasing spin in the spin region 50 {E 1 ^ 1400 keV) is the decay-out effect. Note 
that with the decay-out effect we obtain overall agreements between the calculated results 
and the experimental data. The agreement is satisfactory although it is not perfect. We may 
need more precise comparison using simulated gamma-ray cascades if we want to evaluate 
the differences between the theories and the experiment. 



Finally we show in Fig. 13 the calculated number of bands for 9 Hg. We plot here 



four kinds of result; the calculations where no adjustment of the ND level density nor of the 
collective mass is introduced but we evaluate iVband in two ways (the case (b) as well as the case 
(a)), the calculation where the renormalization of the ND level density with C p = 2.0 x 10~ 4 is 
taken into account (in the case (a)), and the calculation with the adjustment of the collective 
mass (C mass = 3.0) instead of the ND level density. In parallel to the results shown in Fig. |l(], 
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Figure 13: The calculated iVb an d for 192 Hg. The results with C p = C mass = 1 using formula 
(a) and (b) are plotted by a thick solid line and a thin dotted one, respectively. Those using 
(a) with the fitted C p and the fitted C mass are plotted by a thick dashed line and a thick dot- 
dashed one, respectively. They are compared with the number of bands without considering 
the decay-out effect (a thin solid line). 



significant differences in iVband are seen among these calculations. Concerning the calculation 
without any adjustment (C p = C mass = 1), the effect of decay-out is significant and iVb an d 
takes a very small number. However, since this calculation fails to reproduce the decay-out 
of the yrast SD band, it should be an underestimate. On the other hand, when we make 
the adjustments to reproduce the yrast decay-out, the results depend very strongly on the 
method of adjustment. When we adjust the level density of ND states (C p = 2.0 x 10 -4 ), 
the resultant A^and takes the value of more than 100 for / <; 20. This is because the decay- 
out occurs only at very low spins, 1^20, even at high excitation energy (Fig. [[(]), and the 
large value of iVband originating from the shell effect [19] remains. When we perform the 
adjustment with the collective mass (C maS s = 3.0), on the other hand, the decay-out effect 
influences iVband m the whole interval of spin, and reduces iVband significantly taking a value 
less than 70. 

It has been found in Ref. [19| that an extraordinary large iVband due to the shell effect on 
the rotational damping has been predicted (also shown in Fig. 13). However, the present result 
with including the decay-out process suggests the dramatic reduction if we use a realistic ND 
level density and believe the behaviour of the potential energy surface at higher-spin region. 
This is mainly because the barrier height of 192 Hg is lower at higher spins and consequently 
the decay-out boundary is relatively low in excitation energy as is discussed in § 3.2.5| . This 
fact clearly shows that the effect of barrier penetration is crucial in this nucleus, and suggests 
a possibility to study a large-scale change of the potential energy surface as a function of 
angular momentum. Experimental information on iVband would be of crucial importance in 
this respect. 
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4 Conclusions 



We have constructed a microscopic model of thermally excited superdeformed states that 
can describe both the rotational damping caused by the residual two-body interaction and 
the decay-out associated with the barrier penetration to normal deformation. Combining 
the cranked Nilsson-Strutinsky model for the deformed rotating mean-field, the shell-model 
description of individual excited superdeformed states, and the the pair hopping mass for the 
quadrupole shape transition, we define the model microscopically without any phenomeno- 
logical fitting. The calculations have been done for representative superdeformed nuclei 152 Dy 
and 143 Eu in the A ~ 150 region, and 192 Hg in the A ~ 190 region. 

The model thus constructed is applied to the decay-out of the yrast superdeformed states 
to check its validity. The sharpness of the observed decay-out is described well, and a good 
quantitative account of the decay-out spin is achieved for 152 Dy and 143 Eu. However, the 
model fails to describe the decay-out spin in 192 Hg. We examined possible origins of the 
discrepancy in connection with the level density of normal deformed states and the collective 
mass parameter. 

We discussed in detail decay-out properties of the excited superdeformed states having 
the thermal excitation energy up to several MeV above the yrast state. For this purpose, we 
extended the decay-out theory of Vigezzi et al. [11] by using the prescription of Bj0rnholm 
and Lynn 31]. The calculation shows that the excited superdeformed states decay out at 
much higher spins than the decay-out spin of the yrast SD bands. 

In terms of the decay-out and the onset of the rotational damping, the excited superde- 
formed states are classified into three groups having different transitional properties; the 
levels that decay with rotational E2 transitions and form the rotational band structure, the 
levels that decay with rotational E2 but exhibit the rotational damping, and the others that 
mix strongly with ND compound levels and decays to ND levels. The model predicts that 
both the decay-out and the rotational damping influence significantly the effective number of 
superdeformed bands, an observable for the excited superdeformed states obtained recently 
by means of the quasi-continuum spectroscopy. The decay-out causes a characteristic de- 
crease in the effective number of superdeformed band with decreasing spin. The study of 
thermally excited superdeformed bands gives us an opportunity to investigate a large-scale 
shape dynamics as well as the damping of collective rotational motion. 
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